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Abstract 

By using an Ansatz for the pressure (measured in terms of the bag constant) 
of the hadronic gas in equilibrium, we have formulated a new bag-type model 
for the phase transition, i. e. the extended bag model. This allows one to 
take into account the nonperturbative vacuum effects from both sides of the 
equilibrium condition. We have explicitly shown that our model automati- 
cally provides an isentropic equilibrium deconfining phase transition from the 
quark-gluon plasma phase to the hadronic gas phase with a temperature and 
chemical potential dependent bag pressure at the phase boundary. This makes 
it possible to regulate the total entropy and baryon number in both phases, i. 
e. the flow of the specific entropy per baryon across the phase boundary. We 
have also determined the functional form of the temperature and chemical 

potential dependence of the bag pressure outside the phase boundary. 
PACS numbers: 05.70.Fh, 12.38.Aw, 12.38.Lg, 12.38 Mh 
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I. The quark-gluon plasma (QGP) phase is a necessary step in the evolution of the 
Exited Matter from Big Bang to the present days. Apparently, the only way to study this 
phase of the expansion of the Universe is nuclear (heavy ion) collisions at high energies 
which makes it possible "to recreate conditions akin to the first moments of the Early 
Universe, the Big Bang, in the laboratory" [1]. Because of the confinement phenomenon, 
the nonperturbative vacuum structure must play a very important role in the transition 
from QGP to the formation of the hadronic particles (i. e., hadronization) and vice versa. 
As it was underlined in our papers [2] , any correct model of the nonperturbative effects such 
as quark confinement or dynamical chiral symmetry breaking (DCSB) becomes a model of 
the true QCD ground state (i. e., the nonperturbative vacuum) and the other way around. 
Thus the difference between the perturbative (which is always normalizable to zero) and 
the nonperturbative vacua appears to be necessarily nonzero and finite so that it describes 
the above mentioned nonperturbative phenomena at zero temperature. The existence of the 
finite vacuum energy per unit volume - the bag constant [3, 4] - becomes important for a 
realistic calculation of the transition from hadronic to QGP phases and vice versa at nonzero 
temperature as well. 

There are two main approaches to investigate QGP, namely the resummed finite- 
temperature perturbation theory (effective field theory method) [5-7] and the lattice one 
[8]. The former breaks down after the fifth order in the QCD coupling constant g because 
of the severe infrared divergences in the Braaten-Pisarski-Kapusta (BPK) series in powers 
of g m+2n ln n g, but it smoothly incorporates the case of nonzero chemical potential(s). The 
latter is a powerful nonperturbative tool to calculate equations of state for both phases. 
However, up to now, there are no realistic lattice data available for nonzero chemical poten- 
tials) (for problems to introduce it on the lattice see, for example, recent paper [9]). As it 
was emphasized in Ref. [10], at present the phase transition at finite baryon chemical po- 
tential can be only studied within the phenomenological models. The most popular among 
them are, of course, the bag-type models [11], which differ from each other by modelling 
the equation of state of the hadronic phase [12]. The equation of state of the QGP phase is 
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usually approximated by the thermal perturbation theory as an ideal (noninteracting) gas 
consisting of gluons and massless quarks. The generalizations on massive quarks and inter- 
acting QGP are also possible though this is not a simple task because of the nonperturbative 
contributions to the thermodynamic potential, coming from the scales gT and g 2 T, which 
badly affect the above mentioned BPK series. 

Within the framework of the bag-type models, a phase transition between the hadronic 
gas (HG) and the QGP phases is constructed via the Gibbs criteria for a phase equilibrium. 
So the phase transition is necessarily of the first order. This means that the thermodynamic 
quantities of interest are discontinuous across the critical curve. The entropy content of the 
QGP is much higher than in the hadronic phase, so the transition from QGP to HG phase 
is impossible at fixed temperature and chemical potential(s), since in physical processes 
like the heavy ion collisions the entropy cannot decrease. To make the phase transition 
reversible, it is necessary thus to change the values of temperature and chemical potential(s) 
in order to ensue the Gibbs conditions for the phase transition. In the above mentioned 
paper [10] another interesting way was proposed to construct a continuous transition which 
preserves the Gibbs criteria but leads to a temperature and chemical potential dependent 
bag constant. In particular, it is explicitly shown that the specific entropy per baryon, which 
is an important physical observable since it indicates the strangeness production in heavy 
ion collisions, becomes, in this case, continuous across the phase boundary. 

We have recently proposed a new bag-type model, the so-called extended bag model 
[13]. It complements the standard bag model by an Ansatz for the pressure at the phase 
boundary (measured in terms of the bag constant) which allows one to correctly take into 
account the nonperturbative vacuum effects from both sides of the equilibrium condition. 
Our main goal in this paper is to show that our model automatically provides an isentropic 
equilibrium deconfining phase transition with a temperature and chemical potential depen- 
dent bag pressure at the phase boundary. Thus our model makes it feasible to regulate the 
total entropy and baryon number in both phases, i. e. the flow of the specific entropy per 
baryon across the phase boundary. 
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II. For the reader's convinient and for the sake of future application, let us briefly discuss 
the main features of our model in the case of the noninteracting QGP. The QGP state 
equation determines the dependence of the QGP thermodynamical quantities such as the 
energy density e and pressure P on the thermodynamical variables temperature T and quark 
chemical potentials ///. There exist excellent reviews on the physics of the QGP (see, for 
example, Refs. [14-16]), as well as on the phase transitions in it [17]. The QGP pressure (i. 
e., the thermodynamic potential Q, apart from the sign) is [16] 

P=lf SB T' + ^J- + ^4-B, (1) 

where B is the bag constant (see below), while Nf is the number of different quark flavours. 
In what follows we will consider the values Nf = 0,1,2 since the inclusion of the strange 
(s) quark requires a special treatment [18]. The value Nf = describes the case of the pure 
gluon plasma. Note also that the state equation (1) is derived by neglecting quark current 
masses. 

The constant fsB, entering the equation of state (1), 

fsB = fsB(N f ) = f (N f ) = ^ (? + 7 -N f ) (2) 

is the Stefan-Boltzmann (SB) constant, which determines the ideal (noninteracting gluons 
and massless quarks) gas limit. Obviously, for Nf = it equals to the standard SB constant 
of the ideal gluon gas. Note that the temperature dependence of the QGP equation of state 
is dominant over the dependence on quark chemical potentials. 

The energy density e of the noninteracting QGP can be obtained from the thermodynamic 
potential (1) as follows (P = —Q) 

e = 3P + AB, (3) 

so the bag constant determines, in general, deviation from the ideal gas relation between 
pressure and energy density. Let us make a few remarks in advance. Our calculations for 
the noninteracting QGP are not based on the bag model state equation (3) [11, 17]. The 
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constraint, determining the phase transition, will be obtained with the help of the Ansatz 
which is beyond the bag model and it is general (see below, part III and Ref. [13]). We will 
use numerical values for the bag constant which were obtained from a completely different 
source, namely they were calculated in the zero modes enhancement (ZME) model of the 
true QCD vacuum taking the instanton contributions into account as well [2, 13]. 

III. The Gibbs conditions for the phase equilibrium between hadronic gas (HG) and QGP 
phases at T = T c are formulated as follows [11] 

P h = P q = Pc, T h = T q = T c ; 3/1/ = // c , (4) 

where subscripts h, q and c refer to HG, QGP phases and transition (critical or crossover) 
region, respectively. At the same time, the difference between e q — at T = T c can remain 
finite (nonzero) and determines the latent heat (LH), €lh- Let us remind that P q and e q are 
determined by Eqs. (1) and (3), respectively. 

Let us formulate our primary assumption (Ansatz) now. The state equation for the 
hadronic phase (the left hand side of the equilibrium condition (4)) is strongly model de- 
pendent [11, 12]. However, in any model the pressure Ph at any values of temperature T 
and baryonic chemical potential /z, in particular aXT = T c and /i — /i c , can be measured in 
terms of the above mentioned bag constant, i.e. let us put 

P h (T c ^ c ) = —^—B, (5) 

where the so-called parametric functions b h = b h (T c , ji c ) and a h = a h (T c , fi c ) describe the 
details of the HG phase at the phase boundary. Evidently, they may only depend on the set 
of independent dimensionless variables which characterize the HG phase. For example, 

a h = a h (T c , ii c ) = a h (x c , t c , y c , z c ), 

b h = b h (T c , /i c ) = b h (x c , t c , y c , z c ), (6) 



where 



. P* 1 T-> / 7 x 

x c = — , t c = — -, y c = — , z c = fi c Ro, (7) 
T c T c m 



and m denotes the hadron mass while Rq denotes the radius of the nucleon, so it allows one 
to take into account finite size effects due to hard core repulsion between nucleons (extended 
volume corrections) [19]. These variables are independent and all other possible dimension- 
less variables are obtained by combination of these, for example, RqT c = z c /x c , m/T c = 
x clVci /W B 1 ^ = x c /t c ,etc. Also the set of independent dimensionless variables at the 
phase boundary (7) may be extended in order to treat the HG phase in a more sophis- 
ticated way. However, in any case, the parametric functions should be symmetric, i. e. 
a h (T c ,/i c ) = a h (-T c , -/i c ) and b h (T c , fx c ) = b h (-T c , -fx c ). 

From Eqs. (1) and (4-5) at T = T c and /if — /i c /3, one obtains 

M N,W + + ^ = ^B, (8) 

where we introduced a new "physical" (effective) bag constant as follows 

B = (b h + a h + N f )B, (9) 

and it linearly depends on Nf, as it should be, at log-loop level (see our papers [2] and 
Ref. [13]). In connection with our Ansatz (5) a few remarks are in order. The alternative 
parametrization with respect to Nf, namely P/ l (T c ,/i c ) = (b' h /a' h Nf + \)B leads to the 
effective bag constant as B = (b' h + a h Nf + 1)B. The parameteric functions bh and ah, 
as well as b' h and a' h , as was mentioned above, may, in principle, arbitrarily depend on 
dimensionless variables (7). If, for example, a' h vanishes at /jl c — then the linear dependence 
of B on Nf will be spoiled. In other words, the choosen parametrization guarantees the linear 
dependence of B on Nf and the alternative one does not. 

From now on, T c and /i c will be calculated in terms of B and not that of old B, i. e. a 
definite numerical value will be assigned to B. So we consider B as the physical bag constant, 
while B as an unphysical "bare" one. The bag constant is universal and it represents the 
complex nonperturbative structure of the QCD true vacuum. Thus the proposed Ansatz, 
allows one to treat nonperturbative vacuum effects (which are parametrized in terms of B) 
from both sides of the equilibrium condition (4). However, it still remains dependent on 
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the arbitrary parameter ah- In order to eliminate this dependence, let us normalize the 
thermodynamic potential at the phase boundary in Eq. (8) to the standard SB constant 
(2). This yields 

fsB(N f ) = (N f + a h )f SB (N f ) = f SB (0) at N f = 0, (10) 

and one immediately arrives at ah = a/ l (T c , /i c ) = 1. This is our normalization condition and 
it leads to good numerical results for T c and fx c (see Ref. [13]). So the general constraint (8) 
finally becomes uniquely determined, namely (Fig. 1) 

/sbT . + ^ + ^ = _L_ 5 , (U) 

and consequently allows one to investigate the bulk thermodynamic quantites in the vicinity 
of a critical point. Let us emphasize the important observation that the numerical values of 
T c and ji c calculated from the constraint (11) do not depend on the way how one approx- 
imates the equation of state of the hadronic phase. In contrast to the standard bag-type 
models (differed from each other by modelling the hadronic phase [12]), in the extended bag 
model their values depend only on B which incorporates nonperturbative vacuum effects 
from both sides of the Gibbs equilibrium condition (4) as it has been already emphasized 
above. 

The QGP energy density and pressure, however, remain undetermined with our Ansatz 
(5) at this stage. In terms of B they become 

P h (T c ) = PJT C ) = — -B (12) 

1 ' qK } (N f + l)[(N f + l) + b h \ 1 J 

and because of Eq. (3), 

r(T] - r - 3& fe +4(jV / + l) - 

q{ c) ~ C ~ (N f + l)[(N f + l) + b h ] B - {L6) 

As mentioned above, the unknown bh reflects the fact that the hadronic phase state equation 
is strongly model dependent. The unknown bh is the price one pays to determine the above 
mentioned physical quantities with the help of our Ansatz. Precisely for this reason, the bag 



model state equation (3) plays no role in our numerical investigation of the phase transition 
with the constraint (11) from which T c , as well as /i c , can be derived. Concluding this part, 
let us note that our model cetainly requires the coexistence regime between QGP and HG 
phases since e q (T c ), as given by Eq. (13), explicitly depends on bh which describes details of 
the HG phase at T = T c . 

IV. Let us discuss now some general features of the extended bag model. The bare bag 
constant B in Eq. (9) as well as the hadronic gas pressure Ph in Eq. (12) were determined 
at the phase boundary only, namely 

b h (T c , He) + a h {T c , He) + N f 

and 

p (T ]= h(T c ,(i c ) B 

{a h {T c ^ c )+N f )[b h {T c ^ c )+a h {T c ^ c ) + N f Y 1 ) 

respectively. In both equations the normalization condition a^(T c , fx c ) = 1 is to be used. 

Note that in these expressions (in comparison with Eqs. (9) and (12)) we have only restored 

the dependence on T c and /i c , i. e. using relations (6). Let us emphasize that bh in general 

is not constant. Only for a simple hadronic phase (consisting, for example, of only massless 

pion gas) it is constant [13]. The dependence of the QGP pressure P q (1) on T and /i 

remains, obviously, in this case unchanged, since the derivatives of bare B (14) with respect 

to T and /i disappear. The phase transition with the bag pressure determined by Eq. (14) 

was investigated in our previous work [13]. 

One may go outside the phase boundary in our model as well by considering a T and 

\i dependent parametric functions, bh(T,fj) and ah(T,/i), which, at the phase boundary, 

become bh and 1. So extrapolation of the bag pressure (14) outside the phase boundary is 

given by 

B{T ^^b h{ T,,) + a h{ T,,) + N f - (16) 
Hence the HG pressure (15) also becomes dependent on T and \i determining it outside the 
phase boundary as well, i. e. 



P(Tu)- bh( > T >ri i (17) 

h{ ' W (a h (T,ri + Nf)\b h (T,l*) + a h (T,ri + N f Y 1 J 

which, evidently, at the phase boundary becomes Eq. (12) or Eq. (15). The QGP pressure 
P q (1) becomes 

P. = lfMN t )T* + + ^ - ^ ^ - |^ ^ - ^ . (18) 

A simple relation (3) between the energy density and pressure now does not take place. The 
QGP energy density should be deteremined from a general thermodynamic relation, namely 

*= r (fWg) T -''.. < 19 > 

where P g is given by the previous equation. Then one obtains 

B 

e q = 3P, + 4 



b h (T,n) + a h (T,n) + N f 

(b h (T, fi) + a/j(T, fi) + iV/) 2 V <9T d[i> v ' 

where we replaced fif by /i/3. From this equation one recovers the standard relation (3) when 
the bag constant B does not depend on T and /i. In the same way should be determined 
the hadronic energy density 

*= r @, + 0-*- (21) 

where P^ is now given by Eq. (17). One finally obtains 

€h = MT, ,) + a B h (T, ,) + N f f ( T ^T + "|>< T ' A 

Bb h (T,fi)[b h (T,y) + 2a h (T,fj) + 2N f ] , d d, _ 
~ {a h {T^) + N f nb h {T^) + a h {T^) + N f fV dT + ^d^ ah{1 ^ ] ^ W 

Note, Eq. (22) need not be multiplied by the overall factor which takes into account extended 
volume corrections [13, 19] since our parametric functions automatically incorporate them 
(see Ref. [13] and below). Combining this equation with Eq. (20) at the phase boundary, 
for the latent heat £lh = £ c — £h one obtains 
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^=N^-^ + jN^i T W + ^) a " (T ^ (23) 

where, after taking the derivatives, one needs to put T = T c and ji = /i c . In derivation of 
Eq. (23), we have used Ph = P q and e q = e c at the phase boundary T = T c and /i — /i c and 
also the explicit expression (17) which becomes Eq. (15). Thus to leading order the latent 
heat, as given by Eq. (23), does not depend on the details of the hadronic phase which are 
described by the parametric functions bh(T,/i) and a^(T, /i). Note that the second (next-to- 
leading) term is suppressed approximately by one order of magnitude in comparison with 
the first (leading) term (especially for the physically relevant case of the two light quark 
species Nf = 2). The parametric function a,h(T, /i) is a slowly varying function (it is nearly 
constant in the vicinity of the phase transition due to its above mentioned normalization 
condition). 

Using now our values of the effective bag constant B [13], one obtains the numerical 
values of the latent heat, up to leading order, in our model as follows 

e LH (N f = 0) = 3.12 GeV/fm 3 , 
e LH (N f = 1) = 1.82 GeV/fm 3 , 

e LH (N f = 2) = 1.38 GeV/fm 3 . (24) 

At the same time, the QGP energy density e q (20) depends explicitly on the parametric 
functions bh(T, /i) and a/i(T, /i). Since they describe the details of the hadronic phase, so 
this requires the coexistense regime between QGP and HG phases in our model as it was 
already mentioned above. 

Concluding this part of our work, let us remind the reader that the constraint condition 
(11), determining the critical curve, does not depend on how one goes outside the phase 
boundary. Thus the critical temperature T c and the critical chemical potential /i c do not 
depend on the details of the hadronic phase which are decribed by the parametric functions 
bh(T,fi) and a h (T,fi). In our model their numerical values [13] are mainly determined by 
the vacuum effects in both phases. 
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V. It was proposed in Ref. [10] to construct the isentropic equilibrium phase transition 
from QGP to HG by considering a T and /x dependent bag pressure. Due to the above 
general discussion, this can be easily incorporated in our model. The isentropic equilibrium 
phase transition from QGP to a hadronic gas at fixed T and /i via the ratio specific entropy 
per baryon is determined as 

(-) = (-) > (25) 

where s and ub are the entropy and baryon number densities, respectively. For the QGP 
phase they are defined as 

* = (f)„ < 26 > 

and 

»b = -(^). (27) 
In the HG phase these quantities are defined as follows 

and 

- = Q)r W 

In the extended bag model the parametric functions bh(T,fj) and a/ l (T, /i) can only be 
the functions of the dimensionless variables, i. e 

a h (T,fi) = a h (x,t,r,y,z,q), 

bh(T,fi)=b h (x,t,T,y,z,q), (30) 

where 

X= T' t = ^7^' T=7 r' ' z = t iR 0> <! = —■ (31) 

1 1 1 c 171 fi c 
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Evidently, this extrapolates the above introduced set of independent dimensionless vari- 
ables (6-7) outside the phase boundary. These six variables are independent and all 
other possible dimensionless variables are obtained by a combination of these, for exam- 
ple, RqT = z/x, y — x/y, -g^j = x/t, etc. Note that the solutions should be symmetric, 
i. e. bh(T,fi) = bh(—T, —/i) and ah(T,/i) = a,h(— T, —[/) because of the corresponding sym- 
metry in the pressure. Also, the set of independent dimensionless variables (31) may be 
extended, treating the hadronic phase in a more sophisticated way, etc. However, for our 
main purpose in this work (isentropic equilibrium phase transition) this set of variables is 
completely sufficient. 

The HG and QGP pressures are given by Eqs. (17) and (18), respectively. It is convenient 
to introduce the functions defined as 

/ = f(T, ft) = b h (T, fj) + a h (T, fi) + N f = /(*, t, r, y, z, q). (32) 

and 

g = g{T,n)= y T,/ ^ = g(x,t,r,y,z,q). (33) 

a h{-l ,fl>) + 

Then from definitions (26) and (27), one obtains 
and 

N f rr.2 N f S 9 B 

"''T^+S&S-jfrjprf, (35) 

respectively. Here we again replaced jif by /i/3. 

The entropy and baryon number conservation condition (25), on account of relations 
(34-35) and Eq. (17) on account of definitions (32-33), finally becomes 



TV + Bj _ f§-9§ 

where 
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(36) 



so = psB(N f ) + ^/T) 2 , 

no = ^(/VT) + ^(/VT) 3 . (37) 

This relation (36) itself is completely sufficient to guarantee the isentropic equilibrium tran- 
sition (since it is only one relation for two unknown functions) but it is not sufficient to 
determine both functions / and g and hence the functional dependence of bh and on 
T and /i. One needs a second independent relation between them. We suggest the QGP 
fireball condition (s/n B ) q ~ 50 [20] as such. In terms of functions / and g this is 

r3 ^ + B t^ 50. (38) 

Thus one obtains a system of two strongly coupled nonlinear differential equations in partial 
derivatives for / and g parametric functions which, of course, requires a separate consider- 
ation. 

VI. It is instructive, however, to investigate a simplified case which allows us to compare 
our results with those obtained in Refs. [10] and [21]. As mentioned above, in Ref. [10] 
the isentropic equilibrium condition was constructed by extrapolating outside the phase 
boundary the bag pressure (16) only. In terms of our parametric functions / and g (32-33), 
this means that the dependence on the parametric function g in the isentropic equilibrium 
condition (36) is neglected. Then finally the condition (36) becomes 



s o _ ~§rf_ 



(39) 



In terms of the dimensionless variables (31), from Eq. (39) one obtains 

df df ,df ( df df df\ 

xn T- xn t— - x(s + xn ) — = s [y— + z— + q— J , (40) 



where s and n are functions of variable x only and the explicit expressions can be easily 
obtained from relations (37) on account of the above mentioned substitution (31). It is easy 
to check that the general solution of this equation is 

c-c 1 

f(x,t,T,y,z,q) = Ax- c t C2 T C2 - c ' 1 y C:i z C4 q c - C:i - c ' 4 (l + a 1 x 2 + a 2 x 4 ) 4 , (41) 
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where 

a = Nf a = Nf (42) 

1 ef SB (N f y 2 io8tt2 f SB (N f y 1 ' 

Below, for simplicity, all dimensional as well as dimensionless constants, appearing in the 
solution (41), will be always included in the integration "constant" A and the same notation 
will be retained. Constants C's are the constants of separation. In terms of T and ji it 
becomes 

f(T, /i) = b h (T, /i) + a h (T, p) + N f = A^ +C *- C T^T C -^ (l + a lX 2 + a 2 x 4 )^ . (43) 

At the phase boundary x = x c , on account of the constraint (11), it finally becomes 

f(T c , im c ) = b h (T c , iM e ) + N f + l = A^ +C *~ C T^- C \ (44) 

where we used the normalization condition a/ l (T c ,/i c ) = 1. This expression should be finite 
at the end points (T c , 0) and (0,/i c ) of the phase diagram (T c , /x c ) shown in Fig. 1, in order 
to make the hadronic pressure at the phase boundary also finite as it should be. Thus the 
solution (43) becomes 

/(T, /i) = b h {T, /i) + a h (T, /i) + 1 = AT°- Cl (l + a,x 2 + a 2 x 4 ) . (45) 

So the bag pressure (16) at the phase boundary becomes simply B = B / f \T C , /j, c ) , i. e. it 
remains finite but arbitrary (since the "constant" A may depend on T c and /x c ). Nevertheless 
the condition (25) remains continuous across the phase boundary since the ratio in the right 
hand side of the isentropic relation (39) remains, evidently, finite in the limit x — > x c = fi c /T c . 
This is the result of that feature of our model that nonperturbative vacuum effects from both 
sides of the equilibrium condition are taken care of. Apparently, this makes it possible to 
regulate the total entropy and baryon number in both phases, i.e. the flow of the specific 
entropy per baryon across the phase boundary. 

VII. It is easy to see that the asymptotics of the solution (45) outside the phase boundary 
are determined by the separation constant C— C\ and they do not dependend on the variables 
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y, z and q. Indeed, at x — » 0, which means hot HG, at fixed T) and low baryon density 
(/i — > 0), from (45) one obtains that the value of S in this limit is determined by fixed 
temperature and C — C±, i. e. becomes constant, A . In the same way at x — > oo, which 
means cold (T — ► 0), baryon dense matter (at fixed fx), from (45) one obtains that its value 
in this limit now is determined by fixed ii and again C — C±, i. e. tends to constant, A^. 
Then the qualitative behaviour of the bag constant B (16) as a function of x (or T and //) 
is shown in Fig. 2. 

Thus in our model the bag pressure always remains finite and positive (for the simplified 
case this is certainly so and apparently this will be valid for the general case as well) in 
accordance with its definition as the difference between the perturbative (normalizable to 
zero) and the nonperturbative vacua. This definition remains valid also for a T and ii 
dependent bag constant. In the relevant ranges of T and fi, the bag pressure, as obtained 
in Ref. [10], is slowly varying function, i. e. it is nearly constant. This can be seen from 
Figs. 1 and 2 presented in Ref. [21]. This behaviour of the bag pressure is in agreement 
with our model as shown in Fig. 2. In the above mentioned paper [21], this behaviour [10] 
was criticized since it increases with increasing chemical potential at a fixed T and thus 
defies a physical interpretation. However, the finite zeros of the bag constant (especially as 
a function of T) obtained and shown in Figs. 1 and 2 of the above mentioned paper [21] are 
completely unsatisfactory since they strictly contradict the definition of the bag constant. 
As a function of T and /x, it can only tend to zero (in the worst case) but never assume zero, 
i. e. it must always remain positive. These unphysical finite zero points are artifacts of the 
inconsistent approximation used in that paper. Firstly, it is incorrect to directly compare 
results obtained at different scales of the QCD coupling constant, a s = and the running a s . 
The former should be treated at an infinitely large scale while the latter one - at finite, fixed 
scale, so there is no smooth limit to a s = 0. That is why the noninteracting QGP should 
be investigated separately from the interacting one. Secondly, and more importantly, it is 
well known that the first nonperturbative contribution of order a^J 2 to the thermodynamical 
potential is uncomfortably large even for small a s [22]. So the result obtained on account 
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of only the first perturbative correction of order a s will be completely destorded when the 
above mentioned nonperturbative contribution will be incorporated in the calculation of the 
thermodynamic potential. Thus to treat the noninteracting QGP from the very beginning is 
a selfconsistent approximation [10, 13] while it is not like that to treat the interacting QGP up 
to a s order only [21]. In order to achieve that the approximation should be selfconsistent for a 
running a s as well, it is necessary to consider at least the nonperturbative contributions from 
the scale gT. Finally, in Ref. [21], in comparison with Ref. [10] and present investigation, 
there is no phase boundary at all across which entropy per baryon may be continuous on 
account of a T- and /i- dependent bag constant. 

The extension of our model to the case of a running coupling constant (which has been 
already very tentatively treated in Ref. [23]) and the determination of the nontrivial fluc- 
tuations in the interacting QGP will be the subject of the subsequent paper. 

The authors would like to thank J. Zimanyi, K. Toth, Gy. Kluge, T. Biro, P. Levai, and 
T. Csorgo for useful discussions, remarks and support. 
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FIGURES 

FIG. 1. The phase diagram in the plane (T c ,fi c ) measured in inits of MeV. Here the curve is 
shown for the physically relevant case of the two light quarks N f = 2 and the value 5 1/4 = 300 MeV 
was used (see Ref. [13]). 

FIG. 2. The qualitative bahaviour of the bag pressure (16) as a function of x = n/T. The 
upper line corresponds to Aq > A^, while the lower one to - Aq < A^. 
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